Raman Response of Magnetic Excitations in Cuprate Ladders and Planes 
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An unified picture for the Raman response of magnetic excitations in cuprate spin-ladder com- 
pounds is obtained by comparing calculated two-triplon Raman line-shapes with those of the pro- 
totypical compounds SrCu203 (Srl23), Sri4Cu2404i (Srl4), and La6CagCu24C>4i (La6Ca8). The 
theoretical model for the two-leg ladder contains Heisenberg exchange couplings J« and J± plus 
an additional four-spin interaction J cyc . Within this model Srl23 and Srl4 can be described by 
x ~ J||/J ± = 1.5, x cyc := Jcy C /J± = 0.2, jf r123 = 1130 cm -1 and J x r14 = 1080 cm -1 . The 
couplings found for La6Ca8 are x = 1.2, x cyc = 0.2, and J^ 6CaS = 1130 cm -1 . The unexpected 
sharp two-triplon peak in the ladder materials compared to the undoped two-dimensional cuprates 
can be traced back to the anisotropy of the magnetic exchange in rung and leg direction. With the 
results obtained for the isotropic ladder we calculate the Raman line-shape of a two-dimensional 
square lattice using a toy model consisting of a vertical and a horizontal ladder. A direct compari- 
son of these results with Raman experiments for the two-dimensional cuprates R2CUO4 (R=La,Nd), 
Sr2CuC>2Cl2, and YBa2Cu30s+i yields a good agreement for the dominating two-triplon peak. We 
conclude that short range quantum fluctuations are dominating the magnetic Raman response in 
both, ladders and planes. We discuss possible scenarios responsible for the high-energy spectral 
weight of the Raman line-shape, i.e. phonons, the triple-resonance and multi-particle contributions. 

PACS numbers: 75.40.Gb, 75.50.Ee, 75.10.Jm 



I. INTRODUCTION 

Strongly correlated electron systems in low dimensions 
are of fundamental interest due to their fascinating prop- 
erties resulting from strong quantum fluctuations^^. 
Especially in the case of the high-T c cuprate supercon- 
ductors, the role of quantum fluctuations is heavily de- 
bated. Two-magnon Raman scattering has been proven 
to be a powerful tool to study quantum fluctuations in 
the magnetic sector^iSiLSiS. In contrast to the well 
understood magnon dispersion as measured by inelas- 
tic neutron scatteringiSiii*i2iiii4, the quantitative un- 
derstanding of the two-magnon line-shape in the Raman 
response^ and in the optical conduct ivity^iSiiLiS re- 
mains an issue open to debate. 

Interestingly, in the so-called cuprate ladder sys- 
tems like Srl23 or the telephone-number compounds 
(Sr,Ca,La)i4Cu2404i a prominent peak in the mag- 
netic Raman response is observed at the same en- 
ergy of about 3000 cm -1 as in the two-dimensional 
compound a 19 i 20 ' 22 i 23 . In contrast to the gapless long- 
range ordered two-dimensional compound, the quasi one- 
dimensional two-leg ladders are known to be realizations 
of a gapped spin liquid 21 -. Because the elementary ex- 
citations above this groundstate are triplonsSi, we call 
the corresponding Raman response as two-triplon Raman 
scattering. 

On the one hand, one may expect that the Raman 
response is dominated by short-range, high-energy ex- 



citations, suggesting a certain similarity between lad- 
ders and planes, both being built from edge-sharing CU4 
plaquettes. The peak frequencies are in fact at 3000 
cm . On the other hand, the line shape and in par- 
ticular the peak width strongly varies between different 
compounds. In 2D, the peak width is of the order of 
1000 cm" 1 , in La6Ca8 about 500 cm" 1 , in Srl23 and 
SiT4 only (100 - 200) cm -1 . Due to the observation of 
a very sharp two-triplon Raman line in the spin liquid 
Srl4, Gozar et al. have questioned whether the large line 
width in 2D and the related, heavily discussed spectral 
weight above the two-magnon peak can be attributed to 
quantum fluctuations^. 

In the last years theoretical developments in the field 
of quasi one-dimensional systems, namely the quanti- 
tative calculation of spectral densitiea 25 i 26 i 27 i 28 i 29 i 30 i 31 , 
has led to a deeper understanding of magnetic contri- 
butions to the Raman response of undoped cuprate lad- 
ders. Besides the usual Heisenberg exchange terms the 
minimal magnetic model includes four-spin interactions 
which are 4-5 times smaller than the leading Heisen- 
berg couplingsiiii^LS&S 2 ^. The existence and the size 
of the four-spin interactions are consistent with theoret- 
ical derivations of generalized t-J models from one-band 
or three-band Hubbard modela 34 i 35 i 36 i 37 i 38 i 39 i 40 . 

In the present paper we show that the strong variation 
of the line width can be traced back to changes of the 
spatial anisotropy of the exchange constants. The sharp 
Raman line in Srl4 and Srl23 results from x = 1.5, the 
increased line width in La6Ca8 reflects x = 1.2, and the 



isotropic couping x — 1 for the square lattice yields the 
much larger width observed in 2D. In fact, we obtain a 
quantitative description of the dominant Raman peak in 
2D using a toy model which mimics the 2D square lattice 
by the superposition of a vertical and a horizontal ladder. 
We thus conclude that the dominant Raman peak is well 
described by short-range excitations. 

Besides the dominant two-triplon peak, the large 
spectral weight measured at high energies remains an 
open problem for the cuprate ladders and planes. 
We review possible sources of the high-energy spec- 
tral weight which were suggested in the past, 
e.g. quantum fluctuation a 26 i 41 i 42 i 43 i 44 i 45 , the role of 
spin-phonon interactiori 4 * 4 ii* 4 ii 4 £i 4 £i£2*£i an( j the triple 
resonance£i£2i5&5i. In case of the cuprate planes no fi- 
nal conclusion concerning the origin of the high-energy 
weight can be drawn, but in the case of the cuprate lad- 
ders the spin-phonon coupling and the triple resonance 
can be ruled out. 



II. MODEL 

In Raman scattering multi-particle excitations with 
zero change of the total spin can be measured. Starting at 
T = from an S — ground state the singlet excitations 
with combined zero momentum are probed. The Raman 
response in spin ladders has been calculated by first order 
perturbation theory^ and by exact diagonalizatior. 56 . 
In this work, Raman line-shapes are presented obtained 
from continuous unitary transformations (CUT) using 
rung triplons as elementary excitationsSSiS. The results 
are not resolution limited because neither finite size ef- 
fects occur nor an artificial broadening is necessary. 

For zero hole doping, the minimum model for the mag- 
netic properties of the S = 1/2 two-leg ladders is an 
antifcrromagnetic Heisenberg Hamiltonian plus a cyclic 
four-spin exchange term iJ cyd SLSLS 

H = J± ^ Si,iS 2i ,: + J|| ^ S TjJ ;S Ti i + i + H cyc (la) 

i i,T 
H C yc Jcyc ^ -^(l,i),(2,i),(2,i+l),(M+l) 

(lb) 

i 

^(1,1), (1,2), (2,2), (2,1) = -^1234 = (lc) 

(S 1 S 2 )(S 3 S 4 ) + (S 1 S 4 )(S 2 S 3 ) - (S 1 S 3 )(S 2 S 4 ) 

where i denotes the rungs and t £ {1,2} the legs. The 
exchange couplings along the rungs and along the legs 
are denoted by J± and Jm , respectively. The relevant cou- 
plings modeling Srl23 and Srl4S are illustrated in Fig.^ 
There is also another way to include the leading four-spin 
exchange term by cyclic permutationisSL^ 3 . which differs 
in certain two-spin terms from Eq. (|T|l^. Both Hamilto- 
nians are identical except for couplings along the diago- 
nals if J± and Jn are suitably redefined^. 

At T = the Raman response I(u>) is given by the 
retarded resolvent 

I(uj) = --Im(0\OUoj-H + i5)- 1 0\0) . (2) 

\ I I / 
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FIG. 1: Schematic view of a two-leg ladder (notation as in Eq. 
Q). The circles denote the positions of Cu 2+ ions carrying a 
spin 1/2 each. The crystallographic axes are such that x \\ b 
and y \\ a for Srl23 and x || a and y \\ c for Srl4. 

The observables rung (C lcg ) for magnetic light scatter- 
ing in rung-rung (leg-leg) polarization read in leading 
orde*2M2i 

G icg = 4 c s^(S M S M+1 + S 2i4 S 2 , i+1 ) (3a) 

i 

ru„g = ^ng^ SliiSaj . . ( 3b ) 

i 

The factors A^ s and ^4q UI1S depend on the underlying mi- 
croscopic electronic model. It is beyond the scope of the 
present work to compute them. The results will be given 
in units of these factors squared. In this article we will 
only consider non-resonant Raman excitation processes. 
We discuss which laser energy should be used in order to 
investigate the non-resonant regime. 



III. METHOD 

Technically, we employ a CUT to map the Hamilto- 
nian H to an effective Hamiltonian H e f[ which conserves 
the number of rung-triplons, i.e. [Ho,H c g] = where 
H := H\ { j ll=0 . Jcyc=0 y^&. The ground state of H cS 
is the rung-triplon vacuum. For the response function 
I(oj) the observable O is mapped to an effective observ- 
able e ff by the same CUT. The CUT is implemented in 
a perturbative fashion in x = J\\/J± and x cyc = J cy c/ J±- 
The effective Hamiltonian is calculated up to high orders 
(1-triplon terms: 11 th , 2-triplon terms: 10 th order). The 
effective observable O c g is computed to order 10 in the 
2-triplon sector. 

The resulting plain series are represented in terms of 
the variable 1- A SG /(J|| +J ± )£LS where A SG is the one- 
triplon gap. Then standard Pade extrapolants^ yield 
reliable results up to J\\/ = 1 — 1.5 depending on the 
value of J C yc/ J±- Consistency checks were carried out 
by extrapolating the involved quantities before and after 
Fourier transforms. In case of inconclusive extrapolants 
the bare truncated series are used. We will estimate 
the overall accuracy below by comparing with DMRG 
results^. The Raman line shape is finally calculated as 
continued fraction by tridiagonalization of the effective 
two-triplon Hamiltonian. 



3 



Sectors with odd number of triplons are inaccessible 
by Raman scattering due to the invariance of the two 
observables 0*3 and Ojjg llg with respect to reflections 
about the centerline of the ladder—. Thus only excita- 
tions with even number of triplons matter. Therefore 
the leading contributions to the Raman response come 
from the 2-triplon sector. It was shown earlier that the 
two-triplon contribution is the dominant part of the Ra- 
man response at low and intermediate energies^SiSl. The 
role of the four-triplon contribution for the high-energy 
spectral weight will be discussed at the end of this work. 



IV. CUPRATE LADDERS 

In this part we will compare the theoretically obtained 
two-triplon contributions to the experimental line-shapes 
of the cuprate ladders Srl23 and Srl4. The crystals of 
Srl23 have been grown and measured under the same 
conditions as described in Refs.|?3and[ll while the data 
of Srl4 have been provided by Gozar et al?° . The ex- 
perimental Raman line-shape depends strongly on the 
laser energy because resonant contributions are present. 
This becomes apparent in a strong anisotropy between 
the width of the two-triplon peak in leg and rung polar- 
ization for laser energies lu cxc = 2—3 eV. The width of the 
two-triplon peak in leg polarization is much sharper. For 
laser energies w C xc < 2 eV, the strong anisotropy between 
both polarizations vanishes^. It is therefore important 
to figure out which laser energy has to be used for the 
comparison between the non-resonant theory and the ex- 
periment in order to study the magnetic excitations only. 

The first criterion can be gained from the optical con- 
ductivity as for example given in Ref. |2(J the intensity of 
the two-triplon peak develops in the same way as the op- 
tical conductivity. For the non-resonant regime both, en- 
ergy of the incident and scattered light, should be smaller 
than the charge transfer gap (Srl4: At=iok ~ 2.1 
eV). Thus, we have chosen spectra with laser energies 
w cxc = 1.92 eV in the case of Srl4 and w C xc = 1-95 
eV for Srl23. Luckily, the value of the optical conduc- 
tivity is about 100 i7 _1 cm _1 in the sub gap regime at 
u) = cj cxc — E-2T ~ 1.5 eV— (which is 1-2 orders of mag- 
nitude larger than for the 2D cuprates^) yielding a non 
vanishing intensity of the two-triplon peak. Here 2?2T 
denotes the energy of the two-triplon peak. 

The second criterion arises from the polarization de- 
pendence of the two-triplon peak. Depending on the laser 
energy used one can observe a drastic difference in the 
line shape between the two polarizationsi2i2£. While the 
difference in the line shapes is large for w cxc > 2.1 eV, it 
does almost vanish in the case of w eX c < 2.0 e\22. This 
fits very well to the weak polarization dependence of the 
purely magnetic response as described by Eqs. J3J): for 
x cyc = 0.0 the Raman line-shape is identical in the rung 
and the leg polarization. Small deviations x cyc = 0.2 as 
relevant for the description of SiT23 and Srl4 produce 
small deviations with respect to the symmetry between 
the rung and leg polarization of the ladder. These de- 
viations cannot account for the drastic change between 



the two polarizations as observed for u> cxc > 2.1 eV 19,20 . 
We therefore conclude that the spectra lu cxc < 2.0 eV are 
the best choice in order to compare to a purely magnetic, 
non-resonant theory. 

Now we discuss the dependence of the width of the two- 
triplon peak on the parameters x and x cyc . In Fig. [2 the 
full width at half maximum (FWHM) of the two-triplon 
peak is shown. The overall uncertainty shown as error 
bars in Fig. [5] of the extrapolated two-triplon FWHM 
was determined by comparing to DMRG dataS. 

Let us first consider the case x cyc = 0.0. Here the two- 
triplon width should be identical in both polarizations. 
It can be clearly seen that the numerically obtained re- 
sults reflect this property rather well indicating that the 
uncertainties in the extrapolation are small concerning 
the matrix elements. There is a strong dependence of 
the FWHM of the two-triplon peak on the parameter x. 
The peak sharpens significantly when the ratio x of the 
magnetic Heisenberg exchanges increases (4 times from 
x = 1 to x = 1.5). In the case of x cyc ^ 0, the width 
depends on the polarization. In general, the width in 
(xx)-polarization is larger than in (yy)-polarization. For 
fixed x the FWHM changes at maximum by a factor of 
two when varying a; cyc from to 0.2. 

— 1 1 1 1 1 1 1 1 1 - 

- x=1.0(xx) - 

■ x=1.0(yy) 
▲ x=1 .25 (xx) 
a x=1 .25 (yy) 




FIG. 2: Color online The FWHM of the two-triplon peak 
for x := J\\/J± = 1 (squares), x — 1.25 (triangles) and x = 
1.5 (diamonds) as a function of the strength of the four-spin 
interactions x cyc := J cyc /J±. The Orange (Grey) symbols 
denote (xx)-polarization and the blue (black) symbols (yy)- 
polarization. The solid lines are a guide to the eye obtained 
by spline interpolation. 

In Figs. |3I a-d), the experimental Raman response of 
Srl23 and Srl4 is shown for (xx)- and (yy)-polarization 
(red/black and cyan/grey curves). The spectra of Srl23 
were taken in the same way as described in Ref. yjj. The 
data of Srl4 has been made available by Gozar et al. 20 . 
In addition, theoretically obtained two-triplon contribu- 
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tions are displayed (orange/grey and blue/black). Ex- 
perimentally, the width of the two-triplon peak of both 
materials is almost identical (~ 150 cm" 1 ). Only the po- 
sition of the two-triplon peak is different (~ 3140 cm -1 
for Srl23 and ~ 3000 cm -1 for Srl4) which is a result of 
the slightly different Madelung potentials of both com- 
pounds. 

For modeling the Raman response we assume x cyc — 
0.2 for both compounds. This order of magnitude was 
previously obtained for cuprate ladders by inelastic neu- 
tron scatterings^, by infrared absoptionSL 2 ^, Raman 
respons o 19 ' 32 and theoretical works deriving extended 
low-energy Heisenberg models 3 ^*^. In order to account 
for the FWHM and the two-triplon peak position we de- 
termine x — 1.5 and global energy scales J± — 1130 cm -1 
for Srl23 and Jj_ = 1080 cm -1 for Srl4. It was previously 
argued that in Srl4 a charge order of the chain subsys- 
tem modulates the magnetic exchange in the ladders 3 ^. 
This opens a gap in the Raman response which has a 
large effect on the two-triplon peak for the parameters 
x = 1.2 and x cyc = 0.2 which are appropriate for La6Ca8. 
However, the effect is small for larger x-values because 
the induced gap opens well above the two-triplon peak 
at ~ 3600 cm -1 . The set of parameters used above for 
Srl23 and Srl4 describes quantitatively well the Raman 
response as shown in Tab. [I] and Figs. |3fa-d). Especially 
both polarizations for each material can be modeled using 
only one set of parameters Jj_, x and x cyc . The smaller 
FWHM of Srl23 and Srl4 compared to LaeCaSS 3 - 3 ^ can 
be directly explained by their larger a;- values (see Fig- EJ • 
The coupling constants of Srl4 are in good agreement 
with those obtained by IR absorption measurements^. 
Additionally, our set of parameters yields a spin gap of 
290 cm" 1 for Srl23 and 280 cm^for Srl4 using the un- 
derlying one-triplon dispersion. The latter value is con- 
sistent with the spin gap measured by inelastic neutron 
scattering^. 



V. CUPRATE PLANES 

In this section we calculate the Raman response for 
the undoped two-dimensional cuprate compounds using 
a toy model consisting of two uncoupled two-leg ladders 
(see Fig. 0J) . This is motivated by the fact that the build- 
ing blocks of ladders and planes are edge-sharing CU4 
plaquettes. We expect that the Raman response is dom- 
inated by short-range and high-energy excitations yield- 
ing a certain similarity between ladders and planes. In- 
deed, the positions of the two-magnon peak in the 2D 
cuprates and the two-triplon peak in the cuprate ladders 
are found at almost the same frequency ~ 3000 cm -1 , 
but the FWHM of the two-dimensional compounds is a 
factor of 2-6 larger. We have shown in the last section 
that the FWHM of the two-triplon peak in the cuprate 
ladder compounds strongly varies with x. We therefore 
conjecture that the larger FWHM of the two-dimensional 
cuprates originates from the isotropic coupling x = 1. 
There will be of course deviations at small energies re- 
sulting from the differences between a gapped two-leg 




FIG. 3: Color online Comparison of the magnetic Raman 
response of Srl23 (T = 25 K) and Srl4 (T = 5 K) with 
the theoretically obtained two-triplon contribution. The data 
of Srl4 have been provided by Gozar et aim—, (a) The red 
(black) curve denotes the (xx)-polarization (x || 6) of Srl23 
with a laser excitation energy o; CX c =1.95 eV. The orange 
(grey) curve displays the theoretical two-triplon contribution 
with x = 1.5, x cyc = 0.2 and J± = 1130 cm" . (b) The red 
(black) curve denotes the (xx)-polarization (x || a) of Srl4 
with a laser excitation energy lo C xc = 1-92 eV. The orange 
(grey) curve displays the theoretical two-triplon contribution 
with x — 1.5, Xcyc = 0.2 and J± — 1080 cm -1 , (c) (yy)- 
polarization (y || a) for Srl23. Identical parameters as in (a). 
The cyan (grey) curve displays the experimental data and the 
blue (black) curve the theoretical two-triplon contribution, 
(d) (yy)-polarization (y || c) for Srl4. Identical parameters as 
in (b) and the same colors as in (c). 



ladder and the gapless excitations in the two-dimensional 
compounds. Clearly, a magnon description would be the 
proper starting point to treat the long-ranged ordered an- 
tiferromagnetic state. We think, however, that a triplon 
picture which includes the interactions on the quanti- 
tative level can give a good description of the Raman 
response. A similar treatment in terms of gapped quasi- 
particles already led to an improved agreement between 
theory and experiments^^. 

In the following we will show how to deduce the Ai g 
and Big Raman spectra of the square lattice from those of 
the two-leg ladder. Clearly, one should use x = 1 because 
the square lattice is isotropic (J = J|| = JjJ. Starting 
from the Fleury-Loudon operator the observables Blg 
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FIG. 4: Sketch of two uncoupled spin- ladders. Here one lad- 
der is oriented in y-direction and the other in x-direction. We 
approximate the two-dimensional square lattice by the sum 
of these two uncoupled orthogonal ladders. 



(C Alg ) for magnetic light scattering in Bi g (Ai g ) po- 
larization read in leading order for the two-dimensional 
square lattic o^ 2 ^? 



o L 



- 4),Blg ^ 



S^Sj 



(4a) 



. <ij>,x 



o 



Alp 



A), 



Aig 



E 

<ij>,x 



S S i 



<ij>,V 

E 

<ij>,y 



S.S, . (4b) 



Here < ij >,x (< ij >,y) denotes a summation 
over nearest- neighbors in x-direction (y-direction). The 
parameters Ao.eig and Ao.Aig depend on the underly- 
ing microscopic model and are in general not equalS. 
We approximate the two-dimensional square lattice by 
a sum of two uncoupled two-leg ladders, one oriented 
in x-direction, the other in y-direction. The situation 
is sketched in Fig. 0] The summation over both ladder 
orientations will restore the square lattice symmetries. 
Comparing Eq. with Eq. ® one readily deduces the 
following relations 



£> Blg oc (0 lcg - C7 ung ) 
C Alg cx (C log + O rung ) 



(5a) 
(5b) 



between the relevant observables in the two-leg ladder 
and the two-dimensional square lattice. Note that for 
a^cyc = 0, the Raman response in the Ai g polarization 
vanishes due to the property O log |0) = -O rung |0) 2fi . 
The latter point is consistent with earlier treatments of 
the two-dimensional Raman response. But for a finite 
strength of the four-spin interactions x cyci also the Ai g 
polarization is finit o 43 i 46 . 

The theoretical two-triplon contribution of the Bi g 
(panel a) and the Ai g (panel b) polarization is shown 
in Fig. [3] The parameters used are an isotropic cou- 
pling x = 1 and a strength of the four-spin interactions 
■^cyc = 0-0, Xcyc — — 0.1, and Xcyc 0.2. 

The Bi g polarization displays a symmetric two-triplon 
peak which is dominating the Raman response. The four- 
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FIG. 5: Color online Two-triplon Raman response of the 2D 
square lattice for x — 1. (a) Big-polarization for x C y C = 0.0 
(cyan/grey), x cyc — 0.1 (dashed) and x cyc — 0.2 (blue/black), 
(b) Aig-polarization for x cyc = 0.0 (orange/grey), x cyc — 0.1 
(dashed) and x cyc — 0.2 (red/black). Note the different scales 
for the Raman response in Ai g and Bi g spectra. 



spin interactions shift the whole spectrum to lower ener- 
gies and decrease the total intensity. The FWHM of the 
two-triplon peak is approximately given by the average 
width of the two-triplon peaks of an isolated two-leg lad- 
der in rung and in leg polarization. Thus, the width is 
nearly independent of the value of x cyc . 

The Ai g polarization is almost zero for vanishing a; cyc . 
This again reflects the accurate extrapolation of the ma- 
trix elements. For non-zero .r cvc . a finite Ai g contribution 
is realized. The differences in the line-shape between Ai g 
and Big are a pure effect of different matrix elements. 
Compared to the two-triplon peak in the Bi g polariza- 
tion the Ai g polarization displays a two peak structure 
where the second peak is sharper and at higher energies. 

In the following we will compare the theoretical two- 
triplon contribution to the Raman response with low 
temperature experimental data on R 2 Cu04 (w cxc = 2.71 
eV)£, Sr 2 Cu0 2 Cl 2 (oj cxc = 2.73 eV) 5 , and YBa 2 Cu 3 6+5 
(w cxc = 2.71 eV)2> taken from the literature. 

As discussed in section II VI the laser energy used for 
the experiment is a crucial issue. Analogous to the lad- 
ders one should use spectra of cuprate planes measured 
with laser energies below the charge gap for comparing 
to the purely magnetic theoretical response. But it turns 
out that the optical conductivity is rather low (< 10 
il _1 cm _1 ) below the charge gap A which results in a van- 
ishing intensity of the two-magnon peah^. An analogous 
choice of the laser energies below the charge gap as dis- 
cussed for the cuprate ladders is not possible. Therefore, 
we used data measured with laser energies o; oxc = 2.7 
cV> A. At the energy uj cxc = 2.7 eV the optical con- 
ductivity is quite smooth and uj exc ^ A w (1.7 — 2.0) 
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TABLE I: Comparison of the two-triplon peak between experimental data of different cuprate ladder and plane compounds 
and the theoretical results. 

Experiment Theory (CUT) 



Material 


Peak 


FWHM" 


^CXC 


Ref. 


X 


%cyc 


J± 


FWHM 




[cm" 1 ] 


[cm" 1 ] 


[eV] 








[cm- 1 ] 


[cm" 1 ] 


SrCu 2 3 (xx) 


3120 


150-220 


1.96 


this work 


1.5 


0.2 


1130 


180 


SrCu 2 3 (yy) 


3150 


120-180 


1.96 


this work 


1.5 


0.2 


1130 


140 


Sri4Cu24041 (xx) 


3000 


100-160 


1.92 


20 


1.5 


0.2 


1080 


180 


Sri 4 Cu 2 404i (yy) 


3000 


120 b 


1.92 


20 


1.5 


0.2 


1080 


140 


LaeCagCu24 04i (xx) 


3010 


550 6 


2.41 


23.32 


1.2 


0.2 


1130 


580 


La 6 Ca 8 Cu24 04i (yy) 


2950 


350 6 


2.41 


23.32 


1.2 


0.2 


1130 


350 


Sr 2 Cu0 2 Cl 2 


2950 


800-1100 


2.73 




1.0 


0.2 


1190 


1000 


YBa 2 Cu30 6 +5 


2750 


1000-1150 


2.71 


6 


1.0 


0.2 


1110 


940 


La 2 Cu04 


3170 


950-1150 


2.71 


6 


1.0 
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a Lower limit of exp. FWHM: linear background subtracted from data. Upper limit: no background corrections. 

b Exp. FWHM: linear background subtracted from data because background exceeds almost the two-triplon peak heights. 



e\&£2i21. Simultaneously uj cxc — 2.7 eV coincides with 
the triple resonance at tu rcs ~ A + 8 J. The triple reso- 
nance theory predicts two peaks in the Raman response 
at about 2.8J and 4J. The relative intensity of both 
peaks depends on the laser energy. The second peak is 
strongly suppressed at cj rcs . In that sense this laser en- 
ergy can be assigned to be closest to the non-resonant 

In Fig. experimental data and theoretical contribu- 
tions using x — 1 and x cyc = 0.2 are shown for both 
polarizations. Frequencies are measured in units of J. 
We first discuss the Bi g polarization in Fig. EJa)- We 
have chosen the global energy scale J for all experimen- 
tal curves such that the positions of the experimental 
two-magnon and the theoretical two-triplon peaks match. 
This yields quantitatively reasonable values for these 
compounds. In addition, we find quantitative agreement 
between the experimental FWHM and the theoretical 
FWHM of the two-triplon peak. The values of J and the 
FWHM are listed in Tab.|Jfor all compounds. Note that 
the FWHM for x = 1 is larger than for the anisotropic 
case x > 1 as discussed for the ladder compounds. 

Clearly, there are also deviations between theory and 
experiment. As expected, the low-energy spectral weight 
in the theoretical line-shape is larger compared to the 
experimental curves. This is definitely a consequence of 
approximating the two-dimensional square lattice with 
quasi one-dimensional models. There is also spectral 
weight missing at higher energies above the two-triplon 
peak. Possible explanations will be described below. 

The results for the Ai g polarization (shown in Fig. 
HHb)) are explained next. We used the same global energy 
scales J for the experimental curves as determined from 
the Big polarization above. In order to reproduce the 
maximum intensity of the experiment, we multiplied the 
theoretical curve from Fig. IUb) by a factor 5. This im- 
plies that the microscopic parameters ^4o,Big and j4o,Aig 
are anisotropic. A possible reason for this anisotropy 
could be the restriction to the Fleury-Loudon observ- 
able. An extension of this observable to higher orders in 
t/U (four-spin and next-nearest neighbor two-spin terms) 



gives additional contributions to O m& and A1 ^Si. The 
relevance of these contributions has not been analyzed. 

In the experiment a broad hump is measured. We find 
it very promising that the theoretical contribution dis- 
plays the dominant spectral weight just for these energies. 
However, the line-shape can not be resolved completely 
because the dip in the theoretical curve is not observed 
in the experiment. It originates from neglecting the finite 
life-time effects which are already present in the descrip- 
tion of the isolated two-leg ladder 32,70 being the building 
block of our square-lattice toy model. We conclude that 
at least a part of the experimental Ai g polarization origi- 
nates from the finite four-spin interactions. For x cyc = 0, 
there is no purely magnetic contribution to the Raman re- 
sponse for this polarization. A finite Ai g Raman response 
can be regarded as an evidence for the presence of sizable 
four-spin interactions. This follows entirely from symme- 
try arguments and holds true for the full two-dimensional 
model. At higher energies, spectral weight is missing in 
the theoretical contribution in an analogous fashion as in 
the Big polarization. 

VI. HIGH-ENERGY SPECTRAL WEIGHT 

As shown in Sect. HVl and in Sect. the CUT cannot 
account for the missing high-energy spectral weight when 
comparing to the Raman experiments. Also other theo- 
ries proposed previously like calculations based on spin- 
waves^i^, paramagnonsZ&, Jordan- Wigner fermionsS 
and numerical studies^S44iS^& were faced with the 
same problem. Extended theories including (i) multi- 
particle contributions, (ii) spin-phonon coupling, (Hi) 
two-magnon/triplon plus phonon absorption and (iv) 
triple resonance are necessary in order to describe the 
high-energy spectral weight. 

Most of the publications deal with the two-dimensional 
compounds. Here we will try to review these ideas and 
reexamine them in the light of our new results. Especially 
the quantitative results for the cuprate ladders can give 
new insights in this discussion. 
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FIG. 6: Color online Comparison of the two-triplon Raman 
response to the two-magnon Raman line-shape of Nd 2 Cu04 
(u oxc = 2.71 eV, T = 30 K) and Sr 2 Cu0 2 Cl2 (w CX c = 2.73 
eV, T = 5 K). The Raman data of Nd 2 Cu0 4 and Sr 2 Cu0 2 Cl2 
are reproduced from Refs.0and0. The experimental curves 
are smoothed and their zero position is shifted horizontally as 
indicated by the black horizontal lines, (a) Bi g -polarization: 
The blue (black) curve denotes the two-triplon contribution 
with x = 1 and x cyc = 0.2. The global energy scale J is 
chosen such that experimental two-magnon and the theoret- 
ical two-triplon peak merge. This yields J = 1190 cm^ 1 for 
Nd 2 0uO4 (cyan/grey) and Sr 2 Cu0 2 Cl 2 (dashed), (b) Ai g - 
polarization (theory: red/black; experiment: organge/grey): 
Same notations as in (a). Note that the same magnetic ex- 
change couplings J are used. The Ai g -CUT is multiplied by 
a factor 5 in comparison to the curve in Fig. |5] 



(i) multi-particle contributions One open problem is 
the role of multi-particle contributions to the Raman re- 
sponse, i.e. the four-magnon contribution in the case 
of the square lattice and the four-triplon contribution in 
the case of the two-leg ladder. At this stage no quan- 
titative calculations are available. But it is known that 
the multi-triplon spectral weights are sizable for the two- 
leg ladde»2£i2i. The main effect of the four-spin interac- 
tion on the high-energy spectral weight is a small shift 
from the two-triplon to the multi-triplon channels 70 . But 
this shift is not sufficient to account for the high-energy 
spectral weight as observed in experiments. This was 
also found in treatments for the two-dimensional square 
latticei24iM£. However, the complete magnetic infrared 
absorption spectrum (including the high-energy part) 
of La6Ca8 can be described quantitatively by including 
multi-particle contributions^ 7 ^. Here x cyc does not play 
the dominant role for the high-energy spectral weighlsSl. 

It is therefore plausible that these contributions give a 
noticeable effect also on the high-energy Raman response. 
There are also indications that the spectral weight cannot 
be fully explained in this way. For example, the four- 



magnon spectral weight was shown to be negligible for 
the 2D square lattice^. But the magnon-magnon inter- 
action which was not treated in this calculation could 
enhance the high-energy spectral weight. Also quan- 
tum Monte Carlo calculations which include all magnon 
contributions for the two-dimensional Heisenberg model 
seem to explain only a part of the high-energy spectral 
weight 44 . But finite size effects and inaccuracies of the 
analytical continuation can lead to uncertainties in de- 
termining the high-energy spectral weight. 

(ii) spin-phonon coupling The latter observations sug- 
gest that additional degrees of freedom are important. 
It was argued by several authors that the coupling to 
phonons produces a large amount of spectral weight 
above the two-triplon peak^SiiMiLSSiSI. In one approach 
the spin-phonon coupling modulates the magnetic ex- 
change couplings with a Gaussian distribution. Another 
approach introduces a finite spin wave damping induced 
by the spin-phonon coupling. Both scenarios produce 
a significantly broadened and asymmetric two-magnon 
peak as observed in experiments£2*£i. 

Nevertheless, the consistency of a spin-phonon cou- 
pling as suggested above with experiments is not clear. 
The magnitude of this coupling has to be unrealistically 
large in order to describe infrared absorption dataii Ad- 
ditionally, it was pointed out by Freitas and SingbSi 
that the temperature-dependent correlation length and 
the spin dynamics which agree well with purely magnetic 
models does not leave room for such a coupling 82,83 . 

There are no investigations of the role of spin-phonon 
couplings for the case of the cuprate ladder systems. But 
the FWHM of the two-triplon peak can be quantitatively 
understood within a purely magnetic model as shown 
in Sect. IIVI Thus, we conclude that the spin-phonon 
coupling is not strong in the case of the cuprate lad- 
der compounds. Such a coupling leads to a broad two- 
triplon peak in the same way as for the two-dimensional 
case. This is a contradiction when considering the Ra- 
man response and the infrared absorption of cuprate lad- 
ders simultaneously: on the one hand one needs a larger 
anisotropy between leg and rung coupling (larger x) in 
order to sharpen the two-triplon peak in the Raman re- 
sponse again (see Fig. [5} but on the other hand one can- 
not explain the infrared absorption with an substantially 
increased a 27 i 75 . A strong spin-phonon coupling is there- 
fore in contradiction with the results obtained for cuprate 
ladders. This can be also seen as an indication that the 
same holds true for the two-dimensional compounds^. 

(Hi) two-magnon/triplon plus phonon absorption A 
third alternative explaining the high-energy spectral 
weight uses phonons as possible momentum sinks. Here 
a strong spin-phonon coupling is not necessary. The 
idea is based on the work of Lorenzana and Sawatzky 
for infrared absorptioni^^ 4 ^. It is well accepted in 
the case of infrared absorption measurements on cuprate 
laddersS 7 ^ and planesi 7 ^ that the dominant processes 
are magnetic excitations which are assisted by phonons. 
It was realized by Freitas and Singh that similar pro- 
cesses could be important also for the Raman response in 
cuprate planes^i. In an analogous fashion a two-triplon 
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plus (Raman active) phonon process for the Raman re- 
sponse in cuprate ladders could be important. It can 
be used to transfer spectral weight above the two-triplon 
peak leading to an asymmetric line-shape. It is a diffi- 
cult task to determine the relative strength of this process 
compared to the usual two-triplon scattering. 

(iv) triple resonance Additionally, the triple resonance 
was proposed to account for the high-energy spectral 
weight in the two-dimensional compounds^SiSi^. As 
already stated in Sec. ]V\ the experimental spectra of the 
planes are taken in the resonant regime. It is known 
that the triple resonance scenario yields an additional 
peak above the two-magnon peak. Its intensity depends 
significantly on the energy of the incident light in accor- 
dance with experiments^. In principle, the same effect 
is also present in ladder compounds. But for the laser 
energy w cxc =1.92 eV< A considered for Srl23 and Srl4 
the triple resonance condition is not fulfilled. 

Due to the simplified model used for the 2D system no 
conclusion about the high-energy weight can be drawn 
from our results. Because a large spin-phonon coupling 
and the triple resonance can be ruled out for the cuprate 
ladder systems, the observed high-energy spectral weight 
in the cuprate ladder compounds has to be explained 
most probably by the multi-triplon or two-triplon plus 
phonon contributions. 

VII. CONCLUSION 

The first part of this work deals with the theoretical 
understanding of non-resonant magnetic Raman scatter- 
ing on cuprate two-leg ladder compounds, namely Srl23, 
Srl4, and La6Ca8. Therefore we applied a triplon- 
conserving CUT on a microscopic spin-model which in- 
cludes Heisenberg couplings and additional four-spin in- 
teractions. We studied the two-triplon contribution to 
the non-resonant magnetic Raman response. The domi- 
nating feature of the two-triplon contribution is the two- 
triplon peak which has a characteristic FWHM depend- 
ing on the model parameters x and x cyc . We carefully 
chose the experimental data closest to the non-resonant 
regime and compared them with our theory. 

The key observation we found is that the sharpness of 
the two-triplon peak in Srl23 and Srl4 in comparison 
to La6Ca8 can be explained by the stronger anisotropy 
of the magnetic exchange along the rungs and legs of 
the ladder. Indeed, the two-triplon peak width depends 
strongly on the parameter x. Both materials can be mod- 
eled with the parameters x « 1.5 and a; cyc ~ 0.2 but 
different global energy scales Jj_ s» 1130 cm -1 for Srl23 
and J± k 1080 cm -1 for Srl4. The parameters for Srl4 



are in good agreement with infrared absorptionSLS and 
inelastic neutron scattering^ experiments. We conclude 
that the dominating two-triplon peak of the magnetic 
Raman response in cuprate ladders can be consistently 
explained within the microscopic model. The presence 
of a four-spin interaction of the order of 0.2 J± can be 
viewed as a settled issue. 

In the second part of this article we used the results 
found for the two-leg ladder to describe the magnetic Ra- 
man response of the undoped two-dimensional cuprate 
compounds in Bi g and A± s polarization. The contribu- 
tion to the Ai g polarization is only allowed for finite four- 
spin interactions due to symmetry reasons. We use an 
isotropic coupling x — 1 and x cyc = 0.2 for the compar- 
ison with the experimental data. Convincingly, we find 
quantitative agreement for the two-triplon peak position 
and the two-triplon peak width for several compounds. 
Additionally, a sizable spectral weight is found in the 
Ai g polarization consistent with experiments. We con- 
clude that the processes dominating the magnetic Raman 
response are short-ranged. 

The last part deals with the missing high-energy spec- 
tral weight above the dominating two-triplon peak for the 
case of cuprate ladders and planes. We review possible 
sources of this spectral weight like multi-particle contri- 
butions, the role of spin-phonon coupling, a two-triplon 
plus phonon process and the triple resonance to the mag- 
netic Raman response. We deduced from our results that 
the high-energy spectral weight cannot be explained with 
realistic values for the spin-phonon coupling. 

In summary, our calculations lead to an unified un- 
derstanding of the magnetic Raman response in cuprate 
ladder compounds within a purely magnetic model. A 
strong spin-phonon coupling can be excluded for these 
materials. Additionally, we obtained a convincing quan- 
titative description of the dominating two-magnon peak 
in the Raman response of cuprate planes using a toy 
model consisting of two uncoupled two-leg ladders. This 
suggests that the short-ranged triplon excitations might 
be an alternative starting point for the description of the 
two-dimensional cuprate compounds. 
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